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Abstract 

A gradient recovery operator based on projecting the discrete gradient 
onto the standard finite element space is considered. We use an oblique 
projection, where the test and trial spaces are different, and the bases 
of these two spaces form a biorthogonal system. Biorthogonality allows 
efficient computation of the recovery operator. We analyse the approxi- 
mation properties of the gradient recovery operator. 
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1 Introduction 

The gradient reconstruction is a popular technique to develop a reliable a pos- 
teriori error estimators for approximating the solution of partial differential 
equations using adaptive finite element methods [12l Q3J Q] [U O [2] ■ Recently we 
have presented a gradient reconstruction operator based on an oblique projec- 
tion [5]. The oblique projection operator is constructed by using a biorthogonal 
system. In fact, for the linear finite element in simplicial meshes, this approach 
reproduces the so-called gradient reconstruction scheme by the weighted aver- 
aging [H El [2] ■ We proved that the approximation property of the recovered 
gradient for any finite element space is similar to the one obtained by using 
the orthogonal projection with respect to i 2 -norm [HI [2]- In this article, we 
aim at analysing the approximation property of the recovered gradient in one 
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dimension using the oblique projection. This construction is quite useful in ex- 
tending the weighted average gradient recovery of linear finite elements [3] to 
quadrilaterals and hexahedras. 

Let f2 = (a, j3) with a, j3 6 K and a < (3. Let A = {a = x < Xi < ... < x n — f3} 
be a partition of the interval O. We define the interior of the grid, denoted 
int (A), as 

int (A) ={ii£A : 1 < i < n - 1} . 

We also define the set of intervals in the partition A as {ij}™T , where Ii = 
[xi, Xi+i). Two sets A n and B n of indices are also defined as A n = {ieN : < i < n} 
and B n = A n \ {0, n}, respectively. A piecewise linear interpolant of a contin- 
uous function u is written as 4u £ Vh with 

n 
i=0 

where <pi is the standard hat function associated with the point Xi, < i < n. 
We define a discrete space, 

V h = span {0o, - ,(f>n} C H 1 (0) . 

The linear interpolant of u € H 1 ^) is the continuous function defined by I^u — 
Xir=o u ( x i) ( t ) i- However, if we compute the derivative of this interpolant I^u, the 
resulting function will not be continuous. To make the derivative continuous we 
project the derivative of the interpolant, = S"=o u ( x i) ^tt > on t° the discrete 
space Vh- There are two different types of projection. One is an orthogonal 
projection and the other is an oblique projection. The orthogonal projection 
operator, Ph, that projects onto Vh is to find a gh — Ph^r € Vh that 
satisfies: 

3 dx = j n ilt^ dx ' 

Since gh G Vh, we can represent it as an {n + l)-dimensional vector: 

/ 9o \ 

9 = j with g h = y^fft^. 

V 9n ) 



Now the requirement given in equation (1.1) is equivalent to a linear system: 
Mg = /, where M is a mass matrix, and 




with fi = [ —-^dj dx. 
3 Jn dx Vl 
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Here the mass matrix M is tridiagonal. We can reduce computation time greatly 
if we have a diagonal mass matrix. This can be done if we use a suitable oblique 
projection instead of an orthogonal projection. We consider the projection 

Qh ■■ L 2 (n) -► V h , 

which is defined as the problem of finding = Qh^r € Vi such that 

g h X h dx = / < ^ L X h dx, X h E M h , 
Jn dx 

where Mh is another piecewise polynomial space, not orthogonal to Vh, with 
dim (M/j) = dim (Vh), see [8j. In fact, the projection operator Qh is well-defined 
due to the following stability condition. There is a constant f} > independent 
of the mesh-size h such that [SJ [5] 

fr> (J-h4>h dx 

HhWitM < P sup J " „ , 0ft S VJj. 
n h £M h \{0} Wh\\L*{a.) 

In order to achieve that the mass matrix M is diagonal we need to define a new 
set of basis functions for Mh, {Xq, ... , A n }, that are biorthogonal to the standard 
hat basis function we used previously. This biorthogonality relation is defined 
as: 

Xi4>jdx — CjSij, cj 7^ 0, 1 < i, j < n, (1-2) 
where Sij is the Kronecker delta function: 



1, if i = j 
0, otherwise ' 



and Cj is a positive scaling factor. The basis functions for Mh are simply given 
by 



X (x) 
X n (x) 

and for 1 < i < n — 1 



' 2(rr — xi) + (x— xq) 

X —Xl ' 

, otherwise 



%o < x < X\ 



2{x — x n -i) + {x — x n ) 



X n —l ^ X ^ 

otherwise 



■ 2(x—Xi- 1 ) + (x — Xi) 



Xi — ^ *****- X X<i 

i+1 



Xi — + 1 

, otherwise 
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Figure 1: The hat basis function (left) and biorthogonal basis function (right) 
with stepsize h = 1. 



By using an oblique projection Qh the mass matrix will be diagonal. We let the 
diagonal mass matrix be D, so that our system is Dg = /. The values are 
our estimates of the gradient of u at the point x^. So, we estimate the gradient 
by finding g = D^ 1 f, where 

Ixo ^i dx 

We want to calculate the error in this approximation, and find out when gi 
approximates u'(xi) exactly for each %i € A. As in [101 we want to see if gi 
approximates v! (x^) exactly when u is a quadratic polynomial. 



2 Superconvergence 

Theorem 1 Let u £ C°(f2). Then we have 



u(x i+ i) - u(xi-i) . 

-, i € B n , 



X%-\-\ X%—\ 

and 

u(xi)-u{x ) u[x n ) - u(x n -i) 
9o = , 9n = • 

X\ Xq %n &n—l 

Proof: We note that 

SAr) = ^ 2( X -^^H,- Xi) Xi < x < _ foralH e Bn 

, otherwise 
Now, we calculate gi for i € B n : 

Jxq ax h 

Ixo ^ dx ' 
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where 



o,\,,i, = r ( - - ; ) ( 2 (*-**-j dx 



,r 



- - a; l+ i) , 



and 



Therefore, g 



3=0 

= u(xi-i) [ ^ l 1 Xidx + u{xi) ( / -f-X l dx + 

dx \Jxi-l dX 

Jx z dx 

(since <fij and overlap only when j G {i — 1, i, i + 1}) 
_ u(x i+ i) - u{Xi-i) 



x l+1 - n-i 

Now we look at the end-points. We note that 

f Xn ^X dx P™ ^Kdx 

•9° = rk, 7 ^ , » and 9n = r x n J x , ■ 

Jz 90*0 dx J xg (p n X n dx 

Computing as before we get 

u(x\) - u(x ) u{x n ) - u{x n -x) 
9o = , g n = ■ 

X\ Xq X n X n —\ 



We have the following super-convergence in L 2 -norm. This is proved as in [7j GJ] . 



Theorem 2 Let hi = Xi + \ — Xi for < i < n — 1, h — maxo<,<n-i hi, and 
= (xi,Xn-i)- If the point distribution satisfies — hi\ = 0(h 2 ) for 

< i < n — 1. XTien we /icwe i/«e estimate 

<h 2 \\u\\ W :^ { n o) , uew 3oo (n). 

i 2 (O ) 



— - Q ^ Uh 
dx dx 



For the tensor product meshes in two or three dimensions satisfying the above 
mesh condition this theorem has an easy extension. 



2.1 Application to quadratic functions 



Corollary 3 Let u 6 P 2 (R). Then reproduces u' (aij) exactly for all x^ G A, 
w/iere: 



2 ' 
2 



i = 

1 < i < n 

i = n 



and gi = 



* Q = Q x, 



"5 -*fi 



Figure 2: A non-uniform grid with 8 nodes (vertical lines). The points (dots) 
are also shown. 



Proof: We use the result of the previous theorem to get 
_ (axf_ 1 + - axf +1 - bx i+1 ) 

9i — 1 \ 

[Xi-i - x i+1 ) 

_ a ( x l-i ~ x i+i) + b (^i-i ~ x t+i) 

a (Xj-i - Xi+i) {xi-i + Xi+i) + b (xj-i - x l+ i) 
= a (a;j_i + x i+1 ) + b. 



On the other hand, 



u (Xi) = 2a I l+o 

= a (xi-i + Xi+i) + b 
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So, gi reproduces vf (Xi) exactly for i € B n . Now for i = and i = n, we have 
\ (ax 2 + bx\ — ax\ — bx ) a (x\ — Xq) +b(xi — x ) 



9o 



and g n = 



| (El - S ) Xi - Xq 

a (xi - Xo) (Xl + x Q ) + 6 (.'El - Xq) , . 

= a (Xx + Xq) + b, 

Xl - Xq 

\ (axl_ 1 + bXn-i - ax 2 n - bx n ) _ a (xl_ x - x 2 n ) + b (x n _i - x n ) 

o \Xn— 1 ^n) ^ji 
^ ^n) ~t~ >^n) b (x n —\ X n ) . . 

= = a (x„_i + x n ) + b. 

Xfi—l x n 

Since 

u' (Jo) = a (^o + x\) + 6, and v! (x n ) = a (x n -i + x n ) + 6, 
we have go and g n reproduce v! (xq) and v! (i„), respectively, exactly. ■ 



Remark 4 (Uniform Grid) Let A = {a ~ xq < xi < ... < x n — (5} be a uni- 
form grid on the interval ft so that Xi — ajj-i = h, Vi G A n \ {0}, where h is 
some constant, called the stepsize. We note that if our grid is uniform, then 
Xi — Xi ViEi € int(A). So, our gradient recovery operator will reproduce the 
exact gradient of any quadratic function on the interior of a uniform grid. We 
cannot recover the gradients at the endpoints exactly, however, since xq =/= xq 
and x n 7^ x n . 

Corollary 5 Let u £ P2^) with u{x) = ax 2 +bx+c, and let the grid be uniform 
with stepsize h and xq = 0. Then \gi — u'(xi)\ — ah for i = 0, n (i.e. for the 
endpoints of the grid). 

Proof: We will start with the case where i = (i.e. the left endpoint). We 
know from Theorem [3] that go = a (x\ — xq) + b. Since our grid is uniform with 
stepsize h, this simplifies to g = ah + b. u' (xq) = 2ax + 6 = 6, since Xq = 0. 
Therefore 

\g - u (x Q )\ = \(ah + b) - 6| 
= ah. 

The case for i = n (i.e. the right endpoint) is proven similarly. ■ 



For a non-uniform grid, we cannot simplify our approximations using the step- 
size h, since the spacing between each adjacent node is not always equal. We 
did not make any assumption about the uniformity of the grid in Theorem [3j 
Thus \gi — u'(xi)\, i G A n , is not zero for a non-uniform grid This is estimated 
in the following corollary. 
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Corollary 6 Let u (x) = ax 2 + bx + c. Then, 

\gi - u' (xi)\ = a (xi-i + x i+ i - 2a*) Vxi G int (A) , 
1 50 - u' (x )\ = a (xi - x ) , and 
\g n - v! (x n )\ = a(a; n _i - x n ) . 

Remark 7 For i G B n let hi — Xj+i — x^. Then we have 

\9i - u' (Xi)\ = a (xj_x + x i+1 - 2xi) = a(h i+1 - hi). 

We still get superapproximation of the gradient recovery when hi + \ — hi = 0(h 2 ) 
when i G B n . 

2.2 Application to cubic functions 

Corollary 8 Let u G P3QR) with u(x) — ax 3 + bx 2 + cx + d. Then, 
\gi - u' (xi)\ = - (xi-i - x l+ if 

for all i G B n , and 

|.g - u' (i ) I = ^ (x - xxf 

\g n - v! (x n )\ = 2 (x„_i - x n ) 2 

where Xi is defined as in Theorem [3| Similarly, for all i G B n we have 

\g t - u'(xi) \ = a (x$_! + Xi-ix i+ i + x 2 +1 - 3x 2 ) + + x i+ i - 2xi) , 

and 

\g Q - u'(x )\ = a (x\ + x xi - 2x%) + b{x x - x ) , 

\g n - u'(x n )\ = a (a£_i + x n -xx n - 2x 2 n ) + b (x n -i - x n ) . 

Proof: The proof of this theorem is similar to Theorem [3| ■ 

3 Conclusion 

We have presented an analysis of approximation property of the reconstructed 
gradient using an oblique projection. The reconstruction of the gradient is 
numerically efficient due to the use of a biorthogonal system. It is useful to 
investigate the extension to higher order finite elements. 
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